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In this paper we present a model for the kinematics and dynamics of optomechanics 1 which 
describe the coupling between an optical field, here modeled by a massless scalar field, and the 
internal (e.g., determining its reflectivity) and mechanical (e.g., displacement) degrees of freedom of 
a moveable mirror. As opposed to implementing boundary conditions on the field we highlight the 
internal dynamics of the mirror which provides added flexibility to describe a variety of setups rel- 
evant to current experiments. The inclusion of the internal degrees of freedom in this model allows 
for a variety of optical activities of mirrors from those exhibiting broadband reflective properties 
to the cases where reflection is suppressed except for a narrow band centered around the charac- 
teristic frequency associated with the mirror's internal dynamics. After establishing the model and 
the reflective properties of the mirror we show how appropriate parameter choices lead to useful 
optomechanical models such as the well known Barton-Calogeracos model 2 and the important yet 
lesser explored nonlinear models (e.g., Nx coupling) for small photon numbers N, which present 
models based on side-band approximations [3] cannot cope with. As a simple illustrative application 
we consider classical radiation pressure cooling with this model. To expound its theoretical structure 
and physical meanings we connect our model to field-theoretical models using auxiliary fields and 
the ubiquitous Brownian motion model of quantum open systems. Finally we describe the range of 
applications of this model, from a full quantum mechanical treatment of radiation pressure cooling, 
quantum entanglement between macroscopic mirrors, to the backreaction of Hawking radiation on 
black hole evaporation in a moving mirror analog. 
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I. INTRODUCTION 

Optomechanics deals with the interaction of light with mechanical systems. (For an introductory review see e.g., 
P] and references therein.) Though old in name it is relatively new in content - optomechanics has a history at least 
as old as radiation pressure 0J. At the quantum level, optomechanics can be traced at least as far back as Casimir [5], 
who showed that there is an attractive force between two conducting plates from the change of ground state energy in 
the presence of boundary conditions, and to Casimir and Polder 6J who calculated the force on an atom near an ideal 
mirror. The last decade has seen intense interest in several areas that are all under the umbrella of optomechanics. 
To name one such area: the dynamical Casimir effect |7j where a moving object, be it a moving mirror in vacuum, 
a contracting gas bubble in a fluid (e.g., sonoluminescence as advocated in [8]), a time- varying magnetic flux bias 
threading a SQUID terminating a coplanar transmission line [9], or even the spacetime (see below). Optomechanics 
is of renewed current interest because of at least three new developments. The first relating to nanotechnology [10] . 
where miniature mechanical motion can be transduced or manipulated with high precision by capacitive coupling 
or optical control, or in nano-scale wave guides where radiation pressure effects become important, e.g. leading to 
large tailorable photon-phonon couplings which give rise to a large enhancement of stimulated Brillouin scattering 
[llj . The second pertains to quantum information, where information stored in atoms and photons can interface with 
mechanical devices [TH [T3] ■ The third pertains to the use of atoms as optical elements [2] . 

Historically, the gravitation physics community also has explored mirror-field interactions in several ways. For 
example, cosmological particle creation in the early universe (studied by Parker and Zel'dovich in the 60's-70's 
is a form of the dynamical Casimir effect since it arises from the parametric amplification of vacuum fluctuations by 
the expansion of the universe. Another example was the use of a uniformly accelerated mirror as an analog model 
of Hawking-Unruh effects developed by Davies and Fulling [T7J HH]) P30 HOj- It should be noted that these effects 
are different from cosmological particle creation as both the black hole and the uniformly accelerated detector/mirror 
have event horizons while the former in general does not (the de Sitter and anti- de Sitter universes being notable 
exceptions). Yet another example, Forward following the suggestion by Joseph Weber [3T] proposed using laser 
interferometers for the detection of gravitational waves, which has since ushered in today's large-scale and international 
ground-based gravitational wave detection effort [22H24] , These gravitational wave detectors are probably the best 
illustration of the reverse function of optomechanics since in this case an impinging gravitational wave displaces the 
mirrors in the interferometer and the laser beam in the optical arms picks up the corresponding signal. 

In terms of practical applications, an optomechanical process that is actively pursued now is mirror cooling by 
radiation pressure (see e.g., |25]L Optomechanics also provides an excellent means for probing foundational issues in 
quantum physics. Sample studies include: 1) Reaching beyond the standard quantum limit using superconducting 
[26j and nanoelectromechanical |27j devices; 2) Schemes for the improvement of signal-to-noise ratio in gravitational 
wave interferometer detectors [2"BH5T] based on earlier theoretical work of Braginsky and Khalil [3_H i Caves [31 [33J , 
Unruh [34], Kimble, Thorne, et al. [3j; 3) Quantum superposition and entanglement of macroscopic objects such 
as between a mirror and the field [35] and between two mirrors |36[ : and 4) Gravitational decoherence, both in its 
possible limitation to the precision of atom interferometry |37| and as a justification for a modified quantum theory 

Theoretical development for optomechanics also began quite early, most notably in the classic papers of Moore 
[39] . Fulling and Davies [II], Jaekel and Reynaud [40], Barton and Calogeracos [2], Law [41], Dodonov and Klimov 
[412] . Schutzhold, et al. [33] who took a canonical Hamiltonian approach, and Hu and Matacz [33], Golestanian and 
Kardar 45, 46], Wu and Lee [37], Fosco, Lombardo and Mazzitelli 48 who took a path integral approach. There is 
also a lineage of work on relativistically moving mirrors as analog models of the Hawking effect (see e.g., [4"9Tt51| ). 
However, many theoretical aspects remain untouched or were treated loosely (some even erroneously). In view of the 
momentous recent advances in optomechanics, we find it timely and necessary to construct a more solid and complete 
theoretical framework of moving mirrors interacting with a quantum field. 

Our goal is to come up with models and theories capable of treating all of the problems listed above yet conceptually 
simple and theoretically systematic enough to be viable and useful. Admittedly not a simple task [52] , we will delineate 
different aspects as we progress. Suffice it to mention that this first series of papers present the basic models and 
theories of optomechanics both for a closed (this paper) and open (sequel paper in this series) |53| system dynamics 
of moving mirrors in a quantum field. These models can be used to treat the broad class of problems related to the 
dynamical Casimir effect, among other things. The second series includes the back-action of quantum fields on the 
mirror, which is needed for treating mirror cooling (for earlier work see references in [54] ) . the results therein could be 
applied to the related topics of quantum friction [55 -55] and vacuum viscosity [THl 133 US] ■ A third series will focus on 
basic issues in quantum information, making use of the stochastic equations derived in a following paper for moving 
mirrors interacting with a quantum field, specifically on quantum decoherence, superposition, and entanglement of 
mirrors and field. A different vein using similar techniques but staged in curved spacetimes is the moving mirror 
analog of backreaction of Hawking radiation on the evolution of a black hole (for background see references in e.g., 
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EDI). 

In this first paper we tend to the first order of business, namely, that of developing a useful microscopic model 
for any number of mirrors interacting with a field. We consider a massless scalar field in one spatial dimension 
for simplicity. The more realistic electromagnetic field in three spatial dimensions can be treated with a slight 
modification in the form of mirror- field coupling, known as the minimal coupling (see the appendix of |61j). In most 
prior considerations for the primary functions of a mirror its reflective properties (say, by the AMO community) and 
the boundary conditions it imposes on the ambient field (say, by the field theory community), i.e., its amplitude has 
to vanish at the location of the mirror, are considered in a disjoint manner. The advantage of the present model is 
that it avoids the necessity for considering boundary conditions (e.g., a la Fulling and Davies). Only upon elimination 
of the explicit dependence of the internal degrees of freedom of the mirror would the field equations require careful 
attention to boundary conditions. We then consider the kinematics of mirror motion, which also has an effect on 
the field. For example, the motion can parametrically amplify the field modes, including its vacuum fluctuations, 
which results in particle creation (in the field theory language) or "motion- induced" / "acceleration" radiation (in the 
atom-optics language). 

From practical experience physical mirrors have surfaces possessing "light" (as opposed to "heavy") degrees of 
freedom that interact with externally incident radiation in such a way as to maintain the appropriate boundary 
conditions that depend on the material composition of the mirror. Physical mirrors are transparent to sufficiently 
high frequency components of the field because the mirror's internal degrees of freedom are not energetic enough to 
(strongly) couple to field modes with arbitrarily high frequencies. For field modes with frequencies far below this 
cut-off frequency, known as the plasma frequency, the mirror becomes nearly perfectly reflecting. 

In this paper, we treat the mirror motion as that of a particle with mass M and corresponding to the center of 
mass of the mirror. To account for the mirror's reflectivity, we model the mirror's "light" internal degree of freedom 
as a simple harmonic oscillator with mass m and natural frequency 57. This internal variable q(t) is taken to couple 
bilinearly to the massless scalar field at the mirror's location. Because this model involves the mutual interaction 
of the internal oscillator, the field, and the center of mass motion of the mirror we shall call this model an mirror- 
oscillator-field (MOF) model for optomechanical applications. Further details and properties of the MOF model are 
given in Section [TTJ. 

In Section [TT] we demonstrate the mirror's ability to reflect and transmit incident radiation and to perfectly reflect 
or transmit radiation upon judicious choices for the parameter values of the internal oscillator. We also compare 
our model with two commonly used models/descriptions for mirrors: 1) The model of Barton and Calogeracos (BC) 
[2] (described in Section II B I for partially transmitting mirrors; and 2) The auxiliary field approach of Golestanian 
and Kardar [351 HE|- We also show that our model extends the BC model to nonadiabatic regimes of the internal 
oscillator dynamics. In Section |III| we turn our attention to a moving mirror by extending our model to allow for 
arbitrary motion, relativistic or non-relativistic. In Section |IV| we describe the MOF model for multiple moving 
mirrors and focus our attention on how our model appropriately describes multiple reflections and transmissions of 
radiation incident on a cavity. Hence, the MOF model is also applicable to multiple-mirror systems in general and to 
a cavity, in particular, which should be useful for laboratory related studies. In Section |V] we apply the MOF model 
to describe (classical) mirror cooling by radiation pressure and indicate the role of the mirror's internal oscillator. In 
Section |VI| we show how the bilinear coupling in the MOF model relates to the phenomenological model of moving 
mirrors wherein the radiation pressure acts on the mirror through the number of incident photons times the position 
of the mirror (which we refer to as an iVcc-type coupling). In Section VII we show how the MOF model of N moving 
mirrors is related to models of quantum Brownian motion (QBM) involving N harmonic oscillators coupled to a 
bath of oscillators. The available and exact master equations for the latter model will facilitate, among other things, 
our later studies of entanglement between two mirrors, a prototype problem in macroscopic quantum phenomena as 
described above. Finally, in Section [VIII| we summarize our findings and mention further work in progress toward the 
construction of a more complete theory of optomechanics. 



II. A MIRROR AT REST MODELED BY A BILINEAR OSCILLATOR-FIELD COUPLING 



In this section we introduce a model for a mirror at rest interacting with a scalar field. Our system consists of 
a mirror with mass M that we treat as being point-like so that, when allowed to move, its trajectory is described 
by coordinates Z(t). The "light" degrees of freedom, which are responsible for the reflective function of the mirror, 
is modeled as an internal oscillator q(t) with mass m <C M and natural frequency Q. For brevity, we will refer to 
this internal mirror oscillator as a mirosc. Modeling the "light" degrees of freedom by a simple harmonic oscillator 
is functionally similar to the idealization of the internal degrees of freedom of an atom as a "two-level" system when 
considering the atom's optical activities (such as spontaneous and stimulated emission) when interacting with a field 
via a resonant type of coupling |62| from a "harmonic atom" coupling with a bosonic bath when multiple level 
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activities become important. Lastly, we take the mirosc to couple to the external (possibly quantum) scalar field 
<fr(t, x) in a manner that is linear in both quantities (i.e., bilinearly coupled). Taken together, we will refer to this 
model categorically as an mirror-oscillator-field (MOF) model for optomechanical applications. Different oscillator- 
field couplings in this model give rise to different models familiar in optics. However, we will always be considering a 
bilinear coupling in this paper. 

We shall show below how this model can describe, with appropriate choices of parameters, a range of perfectly and 
imperfectly reflecting mirrors. We also show how it relates to the model of Barton and Calogeracos (BC) [2] used in 
the quantum optics community and to the auxiliary field model used more in the field-theory community |45l I46j . 



A. Reflectivity of a mirror modeled by a bilinear oscillator-field coupling 

To demonstrate that the MOF model described above actually possesses the ability to reflect incoming modes it is 
sufficient to put the mirror at rest at the origin so that the action is given by 

S[$,g] = l J d 2 xd a <bd a <5>+ T ^ J dt(q 2 -n 2 q 2 ) + A / tft ?(«)$(*, 0) (2.1) 

where 7y Q( 3 = diag(l,— 1) is the metric of 1+1 dimensional Minkowski space-time. In units where h = c = 1 the 



coupling constant A has dimensions of (mass) 2 = (length) 2 . The equations of motion are obtained by varying (2.1 ) 
in the usual way, 

d a d a $ = <9 2 $ - <9 2 $ = Xq(t)6(x) (2.2) 
mq(t) + mn 2 q{t) = \<f>(t, 0). (2.3) 

Let a plane wave with frequency u) be incident on the mirror from the left (L) so that the field is given by 

$ wi = e ~ tut [9{-x) (e iux + R{uj)e- lulx ) + 9(x)T{u;)e lulx ] (2.4) 

where R(u>) and T(u>) are the frequency-dependent reflection and transmission coefficients, respectively. For the 
steady-state evolution of the mirosc-field system we can take q to oscillate with the same frequency as the incident 
radiation so that 

q(t) = Ae- tut (2.5) 



where the amplitude A is determined from (2.3 1 to be 

ll A — co z to 

The field is continuous at the location of the mirror, ^l^j + ) = $ U L(t, 0~), and the discontinuity of the spatial 



derivative is found by integrating (2.2) over a vanishingly small interval encompassing the mirror's position, 

- dA L (t, 0+) + dA L {t, CP) = Xq(t). (2.7) 



Using these conditions and the equations of motion in (2.2) and (2.3), the reflection and transmission coefficients are 



given in terms of the three mirosc parameters (to, f2, A) by 

R{u) = ~2mu(^-^)+iX» (2 - 8) 
, , 2mw(fi 2 - uj 2 ) . . 

T w = 9 7c^ ^TT^r 2 - 9 

There are three ways that the particle can perfectly reflect incident radiation: 1) In the strong oscillator-field 
coupling limit, A — > oo; 2) When the oscillator is resonantly excited by the monochromatic radiation, 57 = u>, 
(independently of the values of A and to); and 3) In the limit that the mass of the mirosc vanishes, m — > 0. In all 
three cases, the reflection and transmission coefficients are R(lj) = — 1 and T(w) = 0, respectively. 

Likewise, perfect transmission can be attained in three manners: 1) In the limit of vanishingly small oscillator-field 
coupling, A — > 0; 2) When the mirosc frequency is arbitrarily large, fl — > oo; and 3) When the mirosc mass is arbitrarily 
large, m — ^ oo. In all three cases, R(uj) = and T(ui) = 1. 
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FIG. 1. Complex norm of the reflectivity (solid) and transmissivity (dashed) for the case when the minimum in |i?(o;)| 2 is 50% 
(top) and when the oscillator reflects incident modes in a narrow (middle) and a broad bandwidth (bottom). 



The reflection and transmission coefficients possess interesting features that depend on the parameters of the theory 

is 



(m, f2, A). The squared complex norm of the reflection coefficient from (2 



— — I - 7 (2.10) 

—tf-vQ--v > 

where we have introduced the dimcnsionlcss quantity y = oj/Q. To characterize the dependence of the reflection 
coefficient on the frequency oj of the incident field mode we observe that the local minima and maxima occur for 
Umin — l/v3 and y max = {0, 1}, respectively. The maximum at u> = is an artifact of the monopole coupling between 
the field and the mirosc. For a dipole coupling the reflection coefficient vanishes at oj = 0. The reflection coefficient 
equals 1 at both maxima and the value of |i?| 2 at y m i n is 

l#(2/min)| 2 = T-^ (2.H) 

where we define the plasma frequency Q p for the partially transmitting mirror to be 

o3/2\2 

and the index r p = n/n p = 4ml7 3 /(3 3 / 2 A 2 ). 

We can use this minimum in the reflected radiation to indicate when the two maxima of |i?| 2 are sufficiently 
separated and distinguishable. While this is subjective we take |i?(y m i n )| 2 = 1/2 to be our defining requirement, 
which fixes r p — 1. The implication is that if r p 3> 1 then the reflection coefficient is sharply peaked about uj = Q. 
Under this condition, the parameters of the mirosc can be tuned to selectively reflect incident radiation in a narrow 
bandwidth centered on oj — fi, which occurs when the mirosc is resonantly excited (or very nearly so) by the incident 
field. Whenever the mirosc parameters are such that r p -C 1 then the local minimum at y m \ n = l/v3 is close to 1 
and the mirosc reflects modes over a broader frequency bandwidth. Furthermore, the mirror will reflect more than 
half of the incident radiation so long as the frequency of the field is less than ~ £1 to a good approximation. Hence, 
if the mirosc mass is made smaller or the oscillator-field coupling constant A is made larger then more modes will be 
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reflected more strongly by the mirror. Fig. ([!]) shows the basic features of the mirror's scattering properties studied 
in this section and contains plots of the reflection coefficient |i?(y)| 2 for r p equal to, much larger, and much smaller 
than 1. 



B. Relation to common mirror models and approaches 



The bilinearly-coupled oscillator-held model introduced above possesses interesting physical limits that relate to 
two well-known and commonly used mirror models. The first model arises when the mirosc evolves adiabatically 
with the field and gives rise to the model of Barton and Calogeracos (BC) for a partially transmitting mirror. The 
second model arises when the mass of the mirosc becomes arbitrarily small, in which case the mirosc serves the 
role of an auxiliary field that relates to the path integral approaches of [ISJ HSJ [3H] , which describe a quantum field 
interacting with a perfectly reflecting mirror(s). In Section VI we also relate the MOF formulation of optomechanics 



to a commonly used model which describes the effects of radiation pressure by invoking a phenomenological coupling 
between the number of photons impinging the mirror and the mirror's position. 



Barton- Calogeracos model 



The BC model has been used quite often in quantum optics and it is worth summarizing its primary properties 
before showing how it can be derived from our MOF model. Much of BC's attention focuses on quantizing the non- 
relativistic limit of the theory where the mirror velocity is much smaller than c. We do not present their results here 
but refer the reader to the original papers of [2] for further details. 

The action for a mirror at rest in the BC model is, in 1+1 dimensions of space-time, 



Sbc[$] = £ J d 2 x.9 ct $9 Q $-7 J di$ 2 (i,.x = 0) 



(2.13) 



where 7 is related to the plasma frequency of the mirror [2J. Extremizing the action gives the equations of motion 

d a d a 1> = 27<S(af)$(t, 0) (2.14) 
The reflection and transmission of a normal mode of the field incident on the mirror from the left (x < 0) is 



= e 



[9{-x) 



R(w)e 



6»(x)T(w)e i ' 



(2.15) 



where R(ui) and T(u>) are the frequency-dependent reflection and transmission coefficients, reflectively, with the 
property that |i?| 2 + |T| 2 = 1. 

We demand that the field be continuous across the mirror &ujL(t, + ) = &u:L(t,0~) and that its derivative satisfy 



d x $ uL (t, 0+) + d x * uL (t, 0") = 2 7 $ wi (t, 0). 



(2.16) 



This jump condition follows from integrating the field equations across the mirror's position at x — 0. Together with 
the field equation these conditions imply that 



R(u) = 
T{u) = 
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LU + «7 
UJ 

uj + i^i 



(2.17) 
(2.18) 



As the parameter 7 becomes arbitrarily large we see that the reflection becomes perfect and the incoming phase of 
the field changes by 7r radians 



lim R(w) = -1 

7— >oo 



(2.19) 



The ability of the BC model to reproduce perfect and imperfect reflection comes from usin g the quadratic interaction 
<fr 2 (t,0). With this specific coupling to the mirror the jump condition across the origin (2.161 is linear in <!> at the 
mirror, which is vital for obtaining the normal mode in (2.15). 



The MOF model in (2.1) can be related, under appropriate conditions, to the BC model. Observe from (2.3) that 



if q(t) evolves adiabatically with time, 



< 1, 



(2.20) 



7 



then the mirosc follows the time-development of the field at the mirror's position 

A 



lit) 



in 2 



*(t,o) 



Substituting this approximation for the oscillator variable into the scalar field equation (2.2) gives 

A 2 



d a d a >S> 



mfl 2 



<f(ar)$(i,0). 



(2.21) 



(2.22) 



Comparing with (2.14) we recover the model of BC by identifying 7 with the parameters of the mirosc and hence to 



the plasma frequency of the MOF model 



A 2 _ 2 
7_ 2mn 2 ~ 3 3 / 2 p ' 



(2.23) 



Therefore, in the limit that the mirosc changes adiabatically the MOF model yields the BC model. 

An equivalent way of connecting to the BC model is to take the mass of the mirosc to zero, to — > but keep the 
quantity mfl 2 = k constant in this limit, which requires the mirosc natural frequency to approach infinity, fl — > 00. 
In this limit, the mirosc also follows the time-development of the field 



9{t) 



A 



*(t,0). 



(2.24) 



The identification with the BC model then follows the same steps as in the previous paragraph and, in particular, one 
finds that 7 = A 2 /(2k). It is worth pointing out that the massless limit to — > here does not imply that the mirror 
is perfectly reflecting as in the previous section. This is because of the additional requirement that mil 2 = k remain 
constant. In fact, the reflection coefficient (2.8) in this limit becomes 



and the mirror becomes perfectly reflecting when A — > 00. 



2ujk + iX 2 



(2.25) 



Through the identification in (2.23) we may attach heuristic physical interpretations to to, Q (or k) and A. In [? 
], Barton and Calogeracos observe that their model is equivalent to a jellium sheet of zero width, i.e., a surface of 
vanishing thickness having a surface current density generated by the motion of small charge elements with charge 
density n s . If these elements have charge n s e per unit area and mass n s m e per unit area then BC find 



27m s e 2 



7 



(2.26) 



Identifying these microscopic variables to those in our MOF model via (2.23) gives the following relationship 

A 2 



Aim s e 2 

TOp 



(2.27) 



This suggests identifying the mirosc field coupling as a charge per unit area, A — > n s e, and k as a mass per unit area, 
k —> n s m e / '(47r). That is, A can be viewed as a surface charge density and k = mil 2 as a surface mass density. This 
interpretation may be useful for developing a similar MOF model for a mirror in 3+1 dimensions. 



2. Models using auxiliary fields 



The MOF model reduces to another well-known description of mirrors if we take the limit to — > 0. In this limit our 
model describes a perfectly reflecting mirror, as discussed earlier, and the action (2.1) becomes 



lim S[Z,q] = \ 



d 2 xd a <f>d a <f> + A / dtq(t)<f>(t,0) 



(2.28) 



The key point is that the mirosc possesses no dynamics in this limit. Thus, the quantity ip(t) = Xq(t) possesses no 
dynamics of its own and can be regarded as an auxiliary field. 



Oscillator — Field — Mirror model 



m — > 



m —J- 0, mil = constant 



Barton — Calogeracos 
model 




Auxiliary Field 
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2mn 2 



FIG. 2. Relationships among our bilinearly-coupled oscillator-field-mirror model of a partially transmitting mirror, the Barton- 
Calogeracos model of a partially transmitting mirror, and the auxiliary field approach for a perfectly reflecting mirror. 



In the path-integral formulation of the quantum theory, the massless mirosc limit gives rise to the following gener- 
ating functional 63 

Kin Z[J] = J V® J V*P cxp^ J (fx d a $d a $> + i J dtif)(t)$(t,0) + i J d 2 xJ$j (2.29) 

Then, noting that the path integral over ip(t) is just the Fourier representation of the Dirac delta functional, 

J VipexpU J <2t^0)$(t,0) j =£[$(i,0)] (2.30) 

it follows that the generating functional 

\im Q Z[J] = J P$5[$(i,0)] expj^ J d 2 x d a <S>d a <S> + i J d 2 xJ<^ (2.31) 

describes a quantum scalar field constrained to vanish at the location of the mirror (only those field configurations 
that vanish at x — will contribute to the path integral) . The vanishing of the field at the location of the mirror is 
equivalent to the perfect reflection of an incident field [45, 46 . 

Our bilinearly-coupled MOF model (2.1) has successfully reproduced two models describing the interactions of a 
field with a mirror at rest: 1) The partially transmitting BC mirror model when q(t) evolves adiabatically; and 2) 
an auxiliary field approach that enforces the field to vanish at the mirror when the mass of the mirosc is vanishingly 
small. In turn, these two models can be related to each other. Specifically, noting that the delta functional above can 



be approximated by a narrow Gaussian it follows that (2.31) becomes 



lim Z[J] 
m—>0 



fv~,{y*t.«r*-«fM M + i f*. M } (2.32) 



which is increasingly more accurate for larger values of 7. Hence, BC falls out from the generating functional approach 
if we smear the delta functional constraint that enforces the field to vanish on the surface of the mirror. Likewise, using 
the action for the BC model in the generating functional formalism gives the perfect reflection limit when 7 — > 00. 
See Fig. |2]) for the relationships among these theories. 



III. A MOVING MIRROR IN THE MOF MODEL 



As pointed out in the Introduction the physics is quite different in the two cases when the mirror is moving 
relativistically compared to the case when it is moving slowly. The former relates to cosmological particle creation 
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and radiation emitted from black holes or in uniformly accelerated detectors in the Hawking-Unruh effects while the 
latter is closer to accessible laboratory situations such as mirror movements caused by the passing of gravitational 
waves in interferometer detectors and mirror cooling from the field's back-action in the form of radiative pressure 
and quantum friction. The MOF model presented here provides a unified framework for treating both, albeit very 
different, situations. For cases when the mirror motion is prescribed such as coplanar waveguides terminated by a 
SQUID O [64] , or when the mirror possesses non-trivial reflective properties [14] our model can meet the needs of 
current experiments by providing a rich set of reflective properties and a tractable formalism capable of providing 
analytical insight. For systems where the mirror motion is dynamically determined by the mutual interaction of the 
mirror's center of mass, it's internal motion, and the field our model provides a computational ease. This simplification 
results from the fact that boundary conditions are not imposed on the field from the outset but determined by a self- 
consistent elimination of the mirror's internal motion. This facilitates the derivation of equations of motion for the 
mirror's mean position which will be adopted in Sec. [V] to describe classical radiation pressure cooling, and in later 
papers in this series to provide a fully quantum mechanical treatment of mirror cooling and the mirror-analog of the 
black hole back-reaction. 



A. Action formulation 



Allowing the mirror to move requires the addition of an extra term describing its motion along the worldline Z fJ, (X) 
where A is an affine parameter and /i = 0, 1. The physics must re main invariant under any reparameterization of the 
mirror's worldline A — > A(A), which requires modifying the action (2.1) for a static mirror in the following way 



Q 2 q 2 VU a U c 



— M j dX^U a U a + X j dX ^JU a U a q{X)§(Z»{X)) (3.1) 

where an overdot denotes differentiation with respect to the worldline parameter, U^(X) — Z^(X) is the 2-velocity of 
the mirror, dX^/U a U a = dr is the invariant proper time element as measured by an observer on the worldline, and 
indices with Greek letters are raised and lowered with the Minkowski metric ?y Q( 3 = diag(l, —1). The field still couples 
bilincarly to the mirosc via the last term of the action so that the reflective properties studied in the previous section 
are retained by the model. The corresponding Euler-Lagrange equations of motion are easily found to be 



d 2 t S> 



Z"(r)) 



d 2 x <$> = Xq{ T )5 2 {x^ 
mq + mVL 2 q = A$(Z^(r)) 
M eS (r)U» = -H{r){rT ~ U»U")d u $(Z»(r)) 



(3.2) 
(3.3) 
(3.4) 



where we conveniently have chosen to parameterize the worldline by the proper time r a t this point since then 
U a U a = 1 and U a U a — 0, which help simplify the expressions. The quantity M e ff(r) in (3.4) is an effective mass for 
the mirror and is given by 



M ea (r) = M+ \mq 2 {T) 



1 



mnV(r) - Ag(T)*(Z"(r)). 



(3.5) 



Notice that the effective mass has contributions from the rest mass of the mirror (M), the energy of the oscillator 
(mq 2 /2 + mtt 2 q 2 /2), and the interaction energy of the mirror-oscillator-field system (— Xq^(Z)). In other words, the 
effective mass is the rest mass of the mirror plus the total internal energy of the mirosc. 

The structure of (3.2) and (3.4) is reminiscent of a field coupled to a scalar point charge, which here is played by 
the time-dependent mirosc amplitude q{t). In 3 + 1 dimensions, such a system exhibits a radiation reaction force 
on the charge proportional to the third time derivative of the particle's position and exhibits the infamous class of 
unphysical runaway solutions in the absence of any external influences. Below, we show that no such unphysical 



solutions manifest in our MOF model here. To show this, we first solve the field equation in (3.2), which gives 



$(x Q ) = A / dr' G(x a -ZV( T '))q(T') 



(3.6) 

where we ignore the homogeneous solution and where the retarded Green's function in 1+1 spacctimc dimensions is 



G{x a - x' a ) = X -Q(t - t')9(a(x a , x' a )) 



(3.7) 
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where a is half of the squared distance between x a and x' a as measured by the straight line (i.e., a geodesic) connecting 
them, namely, 



a(x ,x ) = -(x -x )(x a -x a ) 



(3.8) 



The derivative of the field evaluated on the worldline is then 

A 
2 



a,$(Z"(r))= [ft,*(s a )] t 



-.Z"{t) 



dr' 9{t - T ')d v a(Z»{T), Z»{t'))8(<j(Z»{t), Z^(r')))g(r') (3.9) 

(3.10) 



+ \ / dr' 8 u (r - r')8(r - T')e(a(Z^T),Z^T / ))q(r') 



From (3.8) it follows that since the mirror's worldline is time-like th en a( Z tJ/ (r), Z^(t')) is always positive except at 
t' — t where it vanishes. Hence, the delta function in the first line of (3.101 receives a contribution only at coincidence, 
when r' — r. Likewise, the 6(t — t') in the second line of (3.10) gives support to the integral at coincidence. 

To evaluate the integrals in ( 3.10 ) we will need to determine the behavior of 6(a) and d v (r — t') when r' ~ r. This 
follows by expanding (3.8) around s = r' — r near zero, giving 

„2 „4 



a(^(r),^(r')) = 



-U a (T)U a ( T ) + 0(s 5 ) 



(3.11) 



2 24 

where we have used the identities U a U a — 1, U a U a — 0, and U a U a = —U a U a , which are valid in the proper time 
parameterization of the worldline. Therefore, writing the delta function in (3.10) as a delta function of s and then 
expanding (3.11) for s near zero gives 



6(a(Z*(T),Z»(T'))) = 



6(s) 



l + -U a U a + 0(s z )). 



In addition, the second integral in (3.10) is proportional to 



dsd„s5(s)0(a)q(T + s) = -q(r) [d v s] s=Q . 



(3.12) 



(3.13) 



The important point to note is that the first integral in (3.10) is potentially divergent. However, we will show now 
that no divergence actually manifests. 

To see this, we observe that (3.8) implies d v a(x a : x' a ) = x v — x' u , which, when evaluated on the worldline and 
expanding around s equal zero, yields 



d v a(Z»(r\ Z»(t')) = -sU v (t) - - U v + 0(s 3 ) 



(3.14) 



Note also that the above equation implies that [9„s] s =o = —U v (t) since from (3.11 ) it follows that d v a = sd„s + 0(s 3 ). 
The integral in (3.10) thus becomes 



dMz^r)) = 



X 



ds 6(-s)( - sU„(t) + 0(, 2 )) (g(r) + 0(e)) *M (l + 0(s 2 )) + \u v (r)q(t). 



Evaluating the integral over s and using 26 (—s) = 1 — sgn(s) we find that 

dMZ»(r)) = \u v (r)q(t), 



(3.15) 



(3.16) 



which is finite. In addition, the derivative of the field above, which is proportional to U v , is contracted with -q^ — U ,1 U V 
in (3.4) to get the force on the mirror, thereby giving zero. Hence, the equation of motion for the mirror's worldline 
from (3.4) is simply 



U»(t) = 



(3.17) 



and the mirror moves inertially. The reason for this trivial dynamics is because the field is not generated by any 
external sources and because we have ignored the initial configuration of the field (i.e., homogeneous solutions to the 
field equation (3.2)). Both of these types of sources will impart a non-trivial dynamics for the mirror's motion. 
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B. Hamiltonian formulation 



Here, we provide a Hamiltonian formulation of the MOF model. To do this, we find it convenient to parameterize 



the worldline by the coordinate time t wherein the action (3.1 ) becomes 

~2 



1 



q,Z] = - dfx d a <S>d a <S> + - dt 



Vi -u 2 



M / dt v 7 ! - U 2 + A / dty/l - U 2 q(t)<S>(t,Z(t)) 



where U (t) — dZ jp and from which the Lagrangian is 



L 



1 



dxd a <Pd a <P 



2 \Vl-U 2 

To derive the Hamiltonian H we first identify the conjugate momenta, 

u(t,x) ,)L 

p(t) 
P(t) 



n 2 q 2 Vl-U 2 ) - My/l -U 2 + Xy/l - U 2 q(t)$(t, Z(t)). 



m(t,x) 

dL mq(t) 



9q(t) y/1 - U 2 (t) 
dL M cS (t)U(t) 



dZ(t) ^l-U 2 {t) 
where the effective mass in terms of the conjugate momenta is 

Meff(t) = M + ^ + l -mn 2 q 2 {t) - \q(t)$(t, Z(t)). 
2m 2 



The Legendre transformation of (3.19) yields the Hamiltonian after some algebra 



H=^jdx(U 2 + (d x $f) + ^/P 2 + M 2 s (t). 



For completeness, we give Hamilton's equations of motion 

$ = n 



Xq(t)M eS (t) 



5(x - Z(t)) 



<1 
P 
Z 

P 



P_ 
m 



Afeff(t) 



^P 2 + M 2 s (t) 
P 

^P 2 + M 2 s (t) 
Xq(t)M cS (t) 



(mfl 2 q- \<$>(t,Z(tj) 



^P 2 + M 2 s (t) 



d x $(t,z(t)), 



(3.18) 



(3.19) 

(3.20) 
(3.21) 
(3.22) 

(3.23) 

(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 

(3.29) 

(3.30) 



which can be shown to be equivalent to the Euler-Lagrange equations in ( 3.2 )-( 3.4 ). As discussed in the previous 
section, an external source will be needed to generate non-trivial forces on the mirror. 

Depending on the application, it may be more convenient to work in a reference frame wherein the interaction 
between the field and the mirror's worldline decouple from each other so that the mirror always remains at rest at the 
origin. A transformation to such a non-inertial frame is advocated in [5| and may be useful for canonically quantizing 
the MOF model. However, we will not pursue this representation here. 
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C. A slowly moving mirror in the MOF model 



Under all laboratory conditions to date the speed of the mirror is small compared to c and justifies developing the 
non-relativistic limit of the mirror-oscillator-field model. For example, it was recently demonstrated that film bulk 
acoustic resonators (FBARs) [55] as large as w 0.5mm can be mechanically oscillated up to 3GHz. The corresponding 
speed of the FBAR (having a modulation depth of 10 -8 ) is only v w 4.4 m/s, which is much smaller than c. Thus, for 
laboratory applications, the non-relativistic limit of the MOF action in (3.1) is entirely appropriate. 



The relativistic Lagrangian (3.19) expanded in powers of Z <C 1 and retaining the lowest order contributions in the 
velocity yields 



L=- I dxd a $d a $> + -mq 2 
2 / 2 



-mtfq 2 + ^MZ 2 - V(Z) + \q{t)<S>(t, Z(t)) 



(3.31) 



where we have dropped the term depending solely on the constant mass of the mirror M and V(Z) describes the 
potential energy of the mirror's motion. The related Hamiltonian follows from a Legendre transform of ( |3.31 1 and is 
found to be 



u = ^ l dx(n 2 + $' 2N| 



p 

2m 



2„2 



-mfl q 



2M 



V(Z)-\q(t)$(t,Z(t)) 



(3.32) 



The equations of motion are easily derived from (3.31 1 or (3.32) so we do not give them here. 



IV. MULTIPLE MOVING MIRRORS IN THE MOF MODEL 



In the previous sections we introduced a model for a mirror whose scattering and reflective properties are described 
by an oscillator, the mirosc, coupled bilinearly to the field. In this section we extend the MOF model to include 
multiple spatially separated partially transmitting mirrors that interact mutually via the field. 

The Lagrangian for N moving mirrors (possibly relativistically) with masses M a (a = 1, . . . , N) can be written as 



i r N / -2 

I -i/*M« + ^(- lf * 



N 



N 



a=l 



a=l 



(4.1) 



and the Euler-Lagrange equations of motion follow straightforwardly and are simply given by Eqs. (3.2 1-( 3.4 ) with all 



mirosc and worldline parameters and variables receiving a subscript a to label the mirror. For completeness and for 
later use, the corresponding Hamiltonian is 



H=\J dx(li 2 + (d x <S>) 2 ) + ^P 2 + M 2 s {t) 



(4.2) 



where the effective mass of the mirror has the same interpretation as before (i.e., mirror rest mass plus total internal 
energy) except now the total internal energy includes the energy of all N mirosc's and their interaction energies with 
the field, 



N 



M eS {t)=M + J2 



^ + \m a n 2 a q 2 a (t) \ a q a (t)<P(t, Z a (t)) 



(4.3) 



In the non-relativistic limit, the Lagrangian and the Hamiltonian are 

N 

L=- f dxd a <S>d a $ + E ( \m«ql ~ \rn a n 2 a q 2 a + hl a Z 2 a - V(Z a ) 

J a=l \ 



X a q a {tMt,Z a (t)) 



N 



H=\J dx(U 2 + (d x <£) 2 ) + E (J^ + \m a nlql + 7^+ V(Z a ) - 



(4.4) 



(4.5) 



In the remainder of this section, we investigate the scattering properties of incident radiation on two mirrors at 
rest. The equations of motion for the two-mirror MOF model are (a = 1, 2) 

2 

d^-d 2 x ^ = J2Kqa(t)6(x) (4.6) 



m a q a + m a Vt 2 q a = \ a ®(t, 0). 



(4.7) 
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Let a monochromatic plane wave of frequency u be incident from the left so that 

$ w z(t,a:) =e~ ib *il> uL (x). 



(4.8) 



The part of the mode ^l(i) can be found using the linearity of the field equation from which the superposition 
principle allows us to write the contributions from multiple reflections and transmissions off of and through both 
mirrors as 



iI^l(x) = 9(-x) 



+9(L - x)6{x) 



Ri + T 1 R 2 T 1 ^(_Ri_R 2 )™ j e~ iu 

n=0 ) 

oo oo 



+6(x - L) 



n=0 

oo 



n=0 



T 1 T 2 ^( J R 1J R 2 )™e^ 



The geometric series can be summed for I-R1-R2I < 1 whereby 

00 

^ I-R1R2 



n=0 



(4.9) 



(4.10) 



To find the reflection and transmission coefficients in terms of the incident frequency uj we assume that the mirosc is 
in a steady-state evolution and oscillates at the same frequency of the radiation so that 



q a (t) = A a e- iut . 
The field is continuous at the locations of each mirror 

^w,(o+) = ^i(o-) 

4> u l(L+)=iJ> uL (L-) 

and the discontinuity of the spatial derivative is to be consistent with the source of the field equation 

-^ L (L + )+^ L {L-) = \ 2 A 2 ; 



(4.11) 

(4.12) 
(4.13) 

(4.14) 
(4.15) 



The mirosc amplitudes A\ , A 2 satisfy the mirosc equations of motion so that we have six equations for the six unknowns 
{R a ,T a , A a } (note the subscript a — 1,2). Thus, the reflection and transmission coefficients are 



1^ 

2mi(jj(VL\ — lo 2 ) — i\\ 

Ti = 1 + Ri 



R-2 



i\\ e 2lu)L 



2m 2 uj(ft 2 — lo 2 ) — iA| 
T 2 = 1 + R 2 e~ 2lulL 

and the amplitude of oscillation for the miroscs are 



An 



l + R 2 
mi {n 2 -uj 2 ) \1 - RiR 2 



A2T9 



m 2 (n 2 -Lo 2 ) \l-R x R 2 



Tie 



iujL 



(4.16) 
(4.17) 
(4.18) 
(4.19) 

(4.20) 
(4.21) 



One can check that the identities |i? a | 2 + \T a \ 2 = 1 are indeed satisfied. The incident field mode (4.9) can then be 
written as 



iPul(x) = 6{-x) 

+9{L-x)9(x) 



Ri + R 2 + 2RiR 2 ^_ iux 



1 — R\R 2 



1 — R\R 2 



R 2 e~ wx ) 



9(x) 



T\T* 



1 — RlR2 



(4.22) 
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When the mirror at x = is perfectly transmitting and the mirror at x = L is perfectly reflecting the field mode is 

^l{x) = 9{L - x) (e iux - e*-^-*)) , (4.23) 

which vanishes as x — > L, as expected. In the complementary case when the mirror at x = is perfectly reflecting 
the field mode incident from the left is 

i/>uL (x) =0(-x) (e iux - e~ iux ) (4.24) 

as also expected. Hence, the MOF model describes the partially reflecting and transmitting properties of two, and 
generally more, mirrors. 



V. CLASSICAL MIRROR COOLING WITH THE MOF MODEL 



In this Section, we show how the MOF model can be used to describe mirror cooling within a completely classical 
context. In a following paper, we discuss quantum effects in mirror cooling using the MOF model [54 . 

The setup is as follows. Consider a cavity formed by two mirrors. We take one of the mirrors to be fixed at the 
origin and perfectly reflecting so that the (classical scalar) field satisfies Dirichlet boundary conditions, &(t, 0) = 0. 
As this fixed and perfectly reflecting mirror will, by assumption, possess no dynamics then we will model the second 
mirror by the MOF model. This second mirror possesses a mirosc internal degree of freedom and will be free to move 
in response to the forces imparted by the field. The motion of this second mirror is assumed to be small relative to 
the size of the cavity, L, and to move on a time-scale much longer than all other time scales in the problem. The 
partial reflectivity of the second mirror allows, for example, a laser field, generated by an external source J ext {x a ), to 
couple to the cavity. 



A. Arbitrary bilinear coupling strength 



The MOF Lagrangian for the system described in the previous paragraph is given by Eq. (3.19) 
1 



L = 



777 A/f 

dx(<S> 2 (x a ) ~ <S>' 2 (x a ) + 2J cxt (x a )<i>(x a )) + -{q 2 {t) - n 2 q 2 (t)) + —{Z 2 {t) - n 2 Z 2 (t)) + \q(tm, L + Z(t)) 

(5.1) 



where we have included an external source J ey ±{x a ) for the field and the second mirror (the dynamical one) has 
coordinates x = L + Z(t) and moves within a harmonic potential with natural frequency Qo. The Euler-Lagrange 
equations for the field, the mirosc, and the coordinates of the movable mirror are 



d 2 $(x a ) - d^(x a ) = J cxt (x a ) + Xq(t)S(x -L- Z{t)) 

q\t) + Q 2 q(t) = -$(t, L + Z{t)) 
m 

Z(t) + n 2 Z(t) = ^q(t) d x $(t, L + Z(t)). 



(5.2) 
(5.3) 

(5.4) 



Our first step will be to solve (5.2) for the field and eliminate its appearance in the remaining equations of motion. 

, so that $ is generated by J ex t and by interactions with the 



Assuming that there is no initial field present [6 
remaining degrees of freedom, then the solution to 



5.2| is 



<J>(x Q ) = J d 2 x / G(x a ;x /a )J cxt (x' a ) + X J dt' G(x a ;t' ,L + Z(t'))q(t') 



(5.5) 



where the retarded Green's function G(x a ;x' a ) for the field subject to Dirichlet boundary conditions at the fixed 
mirror is given by 



G{t,x;t',x') = -9{t-t') 
Note that if x = x' = L > then 



\{t-t') 2 -\{x-x'?) - 



\{t-t') 2 -\{x + x') 2 



(5.6) 



G(t, L;t',L) = - [6(t - t') - 6{t -if- 2L)] . 



(5.7) 



15 



Substituting (5.5 1 into the remaining equations (5.3) and (5.4) gives 



A A 2 
q{t) + n 2 q{t) = -F cxt (t, L + Z{t)) + — I dt' G(t, L + Z(t); t' , L + Z{t'))q{t') 
m m ' 



A 



A 2 



Z(t) + il 2 Z(t) = —q(t)d x F ext (t, L + Z{t)) + — q(t) / dt 1 d x G(t, L + Z(t); t', L + Z(t'))q(t') 



M 



where 



M 



F ext (x a ) = / d 2 x'G(x a ;x' a )J ext (x' a ) 



(5.8) 
(5.9) 



(5.10) 



is the propagated external source for the field. 

Next, we solve for the mirosc variable, q(t). At this point we can take advantage of the assumption that Z(t) <C L 
so that the typical amplitude of the mirror's motion is much smaller than the size of the cavity. This implies we can 
write the solution for the oscillator perturbatively as q = q + qi + • • ■ where q n — 0(Z n ). The equation of motion 
for the leading order mirosc dynamics is 



q (t) + n 2 q (t) - — f dt'G(t,L;t',L)q (t') = -F ext (t,L) 

771 J 777 



The solution to (5.11) is given by (again, ignoring homogeneous solutions) 

/oo 
dt' D{t-t')F eXt {t',L) 
-OO 

where the kernel D(r) is found to be 



D(t) = 



dcu ue~ iujT 

V 2mu}(uj 2 - fl 2 ) + i\ 2 (l - e 2l " L ) 



dui 
2k 



(5.11) 



(5.12) 



(5.13) 



The equation of motion for the first order perturbative correction to the mirosc dynamics is 



A 2 



A 



A 2 



q\{t) + n 2 qi (t)- — I dt'G{t,L;t',L)q 1 {t') = -d x F ext (t,L)Z{t) + — \ dt'[Z{t)d x + Z(t')d x ,]G(t, L;t' ,L)q (t'). 

(5.14) 



The right side of (5.14) simplifies somewhat since ( |5.6[ ) implies that 

d x G(t, L; t', L) = U{t - t' - 2L) = d x ,G{t, L; t' , L) 



and so (5.14) can be written as 
A 2 



A A 2 
dt' G{t, L- 1', L) qi (t') = -d x F cxt (t, L)Z{t) + — (Z(t) + Z{t - 2L))q {t - 2L). 



Thus, the solution to (5.16) is given by 

/°° r A 2 

dt' D(t - t') Xd x F cxt (t', L)Z(t') + —{Z{t') + Z(t' - 2L))q (t' - 2L) 
-oo L * 

Next, we expand the equation of motion for the worldline to leading order in Z{t) to find 

MZ(t) + MtilZ(t) = F[Z(t)} 



(5.15) 



(5.16) 



(5.17) 



(5.18) 



where J r [Z(<)] accounts for the external forces and backreaction from the cavity field and mirosc, and is given by 



A 



A 



F[Z(t)] =\q (t) [d x F ext (t, L) + -q Q (t - 2L)j + Xq (t) ^d 2 F CKt (t, L) - -q (t - 2L)j Z(t) 
+Mo(t) ( - po{t - 2V)Z{t) - ^q (t - 2L)Z(t - 2L) - ^q (t - 2L)Z(t - 2L] 



A 



+Xqi(t) d x F cxt (t, L) + -q (t - 2L) )+ —q a (t) qi {t - 2L) 



A 2 
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We find that the general motion of the mirror as influenced by the cavity field is described by a delay integro-differential 
equation. 

The backreaction terms above will be shown to lead to several effects. First, the driven field will build up in 
amplitude inside the cavity formed by the perfect mirror and the mirror-oscillator. This will lead to a spatially 
varying radiation pressure and a shift in the frequency of the mirror's mechanical motion. Next, depending on the 
equilibrium position of the mirror the cavity field can either accept from or donate energy to the mirror's motion 
arising from retardation effects (see e.g. [T] for a detailed explanation of cooling due to retardation). Finally, non- 
Markovian effects will be present which show how the mirror's motion in the past influences its future movements, 
these effects are accounted for in time-delayed and integral terms. 



B. The weak-coupling limit 



As an example application of these equations we will explore mirror cooling in the weak coupling limit i.e. 
A 2 /(mf2 3 ) <C 1. For many systems of physical interest there exists a large separation between the values of the 
cavity frequency and the oscillation frequency f2o, which allows for a multiple time-scale analysis. In the following 
we will assume that the cavity frequency, the mirosc's frequency fi, and the pump frequency fin are all much larger 
than the frequency of the mirror's mechanical motion f2 - Under these circumstances we may time-average the mir- 
ror's equation of motion in (5.311 over the pump period 2-k/Qd. Since the mirror's mechanical motion is very slow 
compared to this pumping time-scale its trajectory can be safely factored out of any time-averaging integrals so that 

1 f T+t 7(t\ r T+t 

((Z(t)(~-)))=~J dt'Z(t')(...)K-^-J^ dt'(-..) (5.19) 

where ((•)) denotes the time average. 

Let the external source of the field be given by 

J ext (a; a ) = A cos SI D t (5.20) 

so that 

/>oo />oo 

F CKt (x a ) = AJ^ d£ J dx G(t,x;t' ,x')cosQ D t' (5.21) 

where we take the initial time to be at t = and at the end of the calculation take the limit ti — > — oo. Performing 
the spacetime integration gives 

F oxt (t, L) = ae- lUDt + c.c. (5.22) 

where 

.4 



202, 



The first two spatial derivatives of F ext evaluated at x = L are similarly evaluated giving 

8 x F cxt (t, L) = a'e- mDt + c.c. (5.24) 

dlF ext (t, L) = in D a'e- innt + c.c. (5.25) 

where 

-> - lA e in ° L . (5.26) 



2n D 

In addition, we can also derive the explicit form for qo(t) given the expression for the external source 



q (t) = XaD(Q D )e 



(5.27) 



Using these expressions we shall evaluate the time-average of (5.18) over the pump period, T — 2-k/VLd. The 



time-average of the terms independent of q\{t) in (5.19) are easily evaluated. However, the time-average of the terms 
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by their Fourier transforms, 

dv duj 



in (5.18) containing qi(t) requires some elaboration. First, we write qi(t) in (5.171 with Z(t) and D(t — t') replaced 

dt' Z{v)b{Lo)e- wt '- lu{t ~ tr) 



<?i(i) = A 



2?r 2tt 



A 



8 x F ext (t',L) + - (1 + e i2 ^)g (i' - 2L) 



(5.28) 



Evaluating the i' integral and then integrating over u) gives 

/OO 7 f 



A 2 
2 

A 2 
2 



1 + e 



aD(n D )e 



*'* + ^-(1 + e l2uL )a* D* (n D )e- l2nDL 



(5.29) 



In the weak coupling limit, A 2 <C mfl 3 , the effect of the cavity field on the (forced) mirror motion is sufficiently small 
that the mirror will continue to oscillate at a frequency nearly equal to fio- Consequently, we expect Z(y) to be 
sharply peaked for frequencies is ^ CIq. Since L is inversely proportional to the cavity period and v ^ Qq then it 
follows that vL <C 1 and Hd 3> v. Therefore, expanding all terms but Z[y) in the integrand of (5.29) for v near zero 
f? 1 



?i (i) « A<{ e- inot D(n D ) 



a' + \ 2 aD(tt D )e 



2Q,r> L 



e inDt D(-n D ) 



X 2 a*D*(n D )e 



-i2n D L 



Z(t). (5.30) 



Lastly, since the mirror moves on a time scale (~ 1/fio) much longer than the round trip travel time of light in the 
cavity (~ 2L) then fl^L <C 1 and we may expand a ll dela y terms about their instantaneous values so that, for example, 
Z(t — 2L) = Z(t) — 2 LZ(t) + 0((ftoL) 2 ). Using (5.30) one may then easily compute the time-average of the terms 
depending on qi(t) in (5.19). Putting everything together, and remembering to expand the delay terms as discussed 
above, we find that (5.18) becomes 



MZ(t) + T(L)Z(t) + M{nl - An 2 (L))Z(t) = F rll d{L) 



(5.31) 



which is simply the equation for a forced, damped harmonic oscillator with mass M, frequency [fig - Afi^L)] 1 / 2 , and 
damping coefficient T(L). Notice that the latter two quantities depend explicitly on the length of the cavity. After 
time-averaging the explicit form for the radiation pressure is given by 



F rad (L) - \ 2 aD{Vt D ) ( a* + (»•/;•(();,:, ' 



(5.32) 



the frequency shift is given by 



MAn 2 (L) = - in D X 2 aD(n D ) 



3A 2 



a*D*(n D )e 



X 2 b(n D )i a' + X 2 aD(n D )e 



a * + X 2 a*D*{n D )cos2n D L ), 



(5.33) 



and the damping coefficient is given by 



r(L) = —\ab(Q D )\ 2 cos2n D L. (5.34) 

Fig. [3] shows the force on the mirror due to the resulting radiation pressure F rSu d(L) (solid line) and the damping 
constant T(L) (dashed line) as a function of the movable mirror's unperturbed position L from the static mirror at 
the origin. The parameter values chosen for these plots are given in the corresponding figure caption. We observe 
that when the cavity is pumped by an external source, the field energy inside builds up and results in a force from 
radiation pressure F Ta d(L) that varies depending on the size of the unperturbed cavity. 

The gradient of the radiation pressure and the Markov approximation of the integral terms, i.e. those terms 
containing q\(t), leads to a shift in the oscillation frequency of the mirror's center of mass motion. These optical 
spring effects are quantified by the term Af2(L) and changes depending on L. The bottom panel of Fig. shows the 
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FIG. 3. Plots of the force on the mirror from radiation pressure (top left, solid line), the damping constant T(L) (top left, 
dashed line), and the fractional change in the mirror's natural frequency (bottom) as a function of the unperturbed cavity 
length L. The dotted vertical line indicates the value of L = 3.295 ms. The dashed line in the bottom panel shows an 
estimate of where the weak coupling approximation begins to break down. The inset shows the fractional change in frequency 
for L £ [3.287, 3.300] ms. The top right panel shows the evolution of the perturbed trajectory as a function of time with 
L — 3.295ms. The parameter values, in units where c = 1, used to make these plots are as follows: (m,M) — (1,2 000) kg, 
(Q ,Q D ,n) = (1,300tt « 942.5,942)Hz, A = 5000s 2 , and A = 10s" 2 . 



fractional change in the mirror's natural frequency [1 — Af^L)/^ 2 ,] 1 / 2 , which can become imaginary precisely where 
the real part goes to zero in that plot. An important point to note here is that our weak coupling approximation is 
valid when |f2 — AQ(L)| <C £Id- For the values indicated in the figure caption and with L = 3.295ms we see that 
the mirror's modified natural frequency satisfies y/Qg — Af2 2 sa 3.64Hz << tip, which is consistent with the weak 
coupling approximation. It is also important to mention that the mirror's motion can become unstable when the 
mirror is to the right of the resonance, namely, the mirror's spring constant, i.e. K = M(£l 2 — A£! 2 (L)), becomes 
negative as shown in Fig. [3] For the parameter values given in the caption of Fig.|3j the moving mirror's motion is 
damped, the top right panel in Fig.|3j and exemplifies the "cooling" aspect of this classical system to dissipate its 
input energy into the cavity field. 



VI. REDUCTION OF MOF MODEL TO MODELS WITH NX-COUPLING 



In the previous section we used the MOF model in the classical regime to describe the damped motion of one 
mirror of a cavity forced by interactions with an external (laser) field. In the corresponding quantum theory of mirror 
cooling, one usually models the interaction between the mirror and the field by the radiation pressure ~ Nx where 
N is the number operator of quanta (photons) impinging on the mirror's surface and x is the position operator of the 
mirror [1, 4, 33]. We will refer to this interaction as "iVz-coupling." 

The basic motivation for this type of interaction can be easily understood by considering the Hamiltonian for a 
single cavity mode of the form, H cav ~ uj cav (L)aJ a, where a (a') is the annihilation (creation) operator for field quanta 
and w cav (L) is the frequency of a cavity mode of size L. Since the frequency of the cavity modes scales as the inverse 
cavity length w cav (L) ~ 1/L, when we allow the cavity length to vary by a small amount x the frequency is perturbed 
to leading order as uj cav (L + x) ~ w cav (L)(l — x/L + ...). For small cavity length changes the Hamiltonian becomes 

#cav ~ Wcav(£)(l ~ x/L)a^CL 

In this section, we show how the quantum MOF model relates to models with iVx-coupling. In doing so, we highlight 
the assumptions that must be made to connect the two models. We thereby demonstrate that the MOF model should 
be an improvement of the oft used background field approximation for the cavity field [51 I32j In particular, the 
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MOF model should be very useful for studying optomechanical systems having low numbers of cavity photons where 
quantum effects can become quite interesting and important. 

Consider a cavity formed by two mirrors. As in SectionjVJ we take the first mirror (at x = 0) to be fixed for all time 
and perfectly reflecting so that the field satisfies Dirichlet boundary conditions at the origin. We assume the second 
mirror to be partially transmitting and dynamical with small perturbations to its equilibrium position at x = L > 0. 
The second mirror will be described by the MOF model. Recall the Hamiltonian (3.32) for a slowly moving mirror 



H 



dx (n 2 (x a ) + (d x $(x a )) 2 - 2J cxt (x a )<S>(x a )) 



P 2 (t) 
2m 



1 P 2 (t) 

^V(t) + + V(Z(t)) A?(i)*(t, L + Z{t)) 

(6.1) 



where we have included an external source J cx t(x a ) for the field. We shall show, by making a number of assumptions, 
that the interaction component of the above Hamiltonian 



A dx q(t)<Z>(t,L + Z{t)) 



(6.2) 



can be reduced to the TVz-coupling. 

The internal physics of the mirror for many standard radiation pressure cooling calculations is accounted for 
phcnomenologically through the introduction of a cavity quality factor which accounts for the dissipation of field 
energy from within the cavity. In distinction, the MOF model accounts for the detailed information of the mirror's 
internal dynamics. We first will solve for the mirosc to find its effect on the mirror's motion. In this way we trade the 
microscopic information about the mirror for a macroscopic description in terms of the mirror's susceptibility, which 
will establish the link between iVcc-coupling and the MOF model. 

The Heisenberg equations of motion for the field (3.2) and the mirosc ( |3.3[ ) variables are 

(6.3) 
(6.4) 



d^(x a ) - d 2 x ${x a ) = Xq(t)S{x -L- Z(t)) + J cxt (x a ) 



q(t) + n 2 q(t) 



A 



&(t,L + Z(t)). 



We can eliminate t he fi eld's explicit dependence on the mirosc q(t) by solving (6.4) and plugging the solution into 
the wave equation (6.3). In the regime where the mirosc evolves adiabatically so that \q\ <C \Q 2 q\ the approximate 
solution to (6.4) is given by 



K K \l 



(6.5) 



where the second term on the right side is a correction to the leading order, instantaneous solution and is due to the 
fact that the full mirosc dynamics is generally non-Markovian. This can be seen from the general solution of (6.4) 



?(*) = Qh(t) 



A 



dt'g Tdt (t;t')${t',L + Z(l!)) 



where qh(t) is the homogeneous solution and <7 re t is the retarded Green's function for the mirosc 

„. .. sinf2(£ — t') 
g lct (t;t') = e(t-t>) K - >-. 



(6.6) 



(6.7) 



More specifically, the mirosc receives contributions from the past as implied by the integral in (6.6). However, the 
approximation (6.5) is valid if the mirosc degree of freedom is "light" thus responding nearly instantaneously to 
external influences. 



Substituting the approximate mirosc solution (6.5) into the wave equation (6.3) gives the effective dynamics for the 
field 



A 



A 2 



d 2 ^{x a ) - d 2 x <S>{x a ) = J cxt (x a ) + — $(£, L + Z(t))5{x -L- Z(t)) - L + Z{t))8(x -L- Z{t)) + 



K 



kQ 2 



(6.8) 



Notice that (6.8) can be derived from the following effective Lagrangian 

v2 



L ^ = \[ dx(d a ^d a ^ + 2J cxt ^) + ^ J dt<S> 2 {t.L + Z(t)) + ^ 2 J dt^ 2 (t,L + Z(t)) + --- + hlZ 2 



V(Z) 



(6.9) 
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where • • • denotes the higher order terms in (6.5). The interaction Hamiltonian corresponding to the above effective 
Lagrangian is found to be 

H cSint = -^ 2 (t 1 L + Z(t)) + - X ^U 2 (t,L + Z(t)) + -- - . (6.10) 



Assuming that Z(t) <C L, we may expand the effective interaction Hamiltonian in (6.10) to find 

H cSint = H^ int + ffl + 0(Z 2 ) (6.11) 

where 

H§\ nt = -^(t'L) + n 2 (*,i) + • • • (6-12) 
H^ inb = - — Z(t)$(t,L)d x $(t,L) + ^Z(t)U(t,L)d x U(t,L) + -- . (6.13) 

Notice that the leading order interaction Hamiltonian is independent of Z(t) so that it exerts no force on the movable 
mirror. In fact, one can group H^ int with the free Hamiltonian for the field that, when taken together, describes the 
free evolution of the field in a cavity where one mirror is fixed at x = and perfectly reflecting and the other mirror 
is fixed at x — L but partially transmitting. The remaining terms in the effective interaction Hamiltonian describe 
the perturbative response of the second mirror to its coupling with the field and vice versa. 

To leading order in Z(t), we can express the field in terms of a homogeneous solution via the cavity's normal modes 
and in terms of the external source J ox t, 

$(x Q ) w Y,M a kMx)e-*" kt + H.c.) + f d 2 x'G c ^(x a ,x ,a )J cxt (x' a ) (6.14) 
k J 

where H.c. is the Hermitian conjugate of the preceding terms, u k (x) are the normal modes of the cavity and satisfy 

A 2 \ 

\+u% + — 5(x- L))u k (x) =0 (6.15) 
such that Ufc(0) = since the mirror at x = is perfectly reflecting. The retarded Green's function here satisfies 



(d 2 x +loI + ^S(x - L)) G^(u k ; x, x') = -5{x - x') 



(6.16) 



with Dirichlet boundary conditions at the origin G^(uj k ; x, 0) = and G^(oJk'i 0, x') = 0. 

Also, Nk is chosen so that x), H(t, x')] = ihS(x — x') for x and x' greater than zero. These commutation 
relations require a k and a\ to be annihilation and creation operators, respectively. For the following we focus entirely 
on the component of the interaction Hamiltonian coming from the field inside the cavity. The field outside of the 
cavity gives rise to a constant and position independent radiation pressure that only yields a shift in the equilibrium 
position of the mirror at x = L. 

If the cavity is pumped by a laser beam with a frequency slightly detuned from one of the cavity resonances and 
if the cavity quality factor is large then the cavity field, represented as a mode sum, can be approximated well by a 
single mode. Expressing the field in terms of the fundamental cavity resonance we find, at linear order Z(t), that the 
interaction Hamiltonian (6.131 is given by 



H 



(i) 



A 5 



Z(t) 



iV fc (a feUfe (L) e - 4Wfct +H.c.) +F cxt (t,L) N k (a k u k (L)e~ lulkt + H.c.) + d x F cxt {t, L) 



(6.17) 



where • • • refers to corrections arising from time derivatives of the field appearing in (6.13) and 



F ext (x a )= / d'x'G c ™(x a -x' a )J ext {x' a ). 



(6.18) 



For many systems of interest the frequency of the fundamental cavity mode is much larger than the typical frequency 
of the mirror's motion (i.e. Qo/ujk <C 1). Under such conditions the mirror's position changes adiabatically over many 
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oscillations of the cavity field allowing a time average (denoted by double angled brackets) of the effective interaction 
Hamiltonian 

((*$)) = m I dtH^. (6.19) 

Here, T is the period of the cavity's fundamental mode and N is a large integer such that (27r)/fio ^ NT. Since Z(t) 
is approximately constant over the entire integration range it can be taken outside of the time-average giving 

((H$)) « -^Z(t)F cxt (t,L)d x F cxt (t,L) - (^-\N k \ 2 u' k (L)uUL)^Z(t)ala k + R.c. (6.20) 

This step is equivalent to taking the rotating wave approximation. The key point is that the first term on the right 
side is a classical radiation pressure originating solely from the external source while the second term is a quantum 
mechanical radiation pressure and is, in fact, the ./Vi-coupling. 

Before concluding this section, we collect the main assumptions used in relating the MOF model to the phenomeno- 
logical radiation pressure interaction Hamiltonian. The assumptions are as follows: 

• The movable mirror is only ever slightly perturbed from its otherwise equilibrium position at x = L: 

• The cavity frequency is much less than the natural frequency of the mirosc; 

• The cavity has a high quality factor; 

• The cavity is pumped by a laser at a frequency slightly detuned from one of the cavity resonances; and 

• The cavity frequency is much greater than the typical timescale associated with the mirror's motion (i.e., the 
natural period if in a harmonic trap) 

Under these assumptions we have shown that the effective interaction between the mirror and the cavity field is given 
by an ./Vx-coupling. It is possible that the iVx-coupling can be obtained using a different setup and assumptions. 
However, our purpose here is not to elucidate all the ways that the iVrr-coupling can be derived from the MOF model 
but rather to show that it can be derived from a microphysics model of a moving mirror. 



VII. MIRROR-OSCILLATOR-FIELD (MOF) MODEL AND QUANTUM BROWNIAN MOTION 

In this Section we shall establish a connection between the MOF model for N moving mirrors and N harmonic 
oscillators interacting with a bath of harmonic oscillators that constitute an environment for the N oscillators. The 
latter system has a long and well-developed history for providing a simple model with which to study quantum 
Brownian motion (QBM). Hence, if a relationship between the MOF model and QBM exists then one should be 
able to exploit the results of many previous studies (regarding decoherence, (dis)entanglement, fluctuation-dissipation 
relations, etc.) to apply towards moving mirror systems. We show here that such a relationship does indeed exist. 



A. Static mirrors and QBM 

Consider a mirror at rest that is fixed at Z(t) = for all time. The MOF Hamiltonian for this configuration follows 



from (3.32) 



H=- dx(ll 2 (x a ) + (d x $(x a )) 2 ) + ^Ai + -mn 2 q 2 (t) - \q(t)$(t, 0). (7.1) 



It is well known that a field can be represented as a continuum of harmonic oscillators, some of which have arbitrarily 
large natural frequencies. However, such large frequencies are not usually physically relevant (and often lead to 
divergences that must be properly handled with well-established renormalization techniques) so that one can simply 
impose a cut-off frequency A, which has the effect of ensuring that all calculated quantities are finite |67j . 
The mode decomposition of the field is 

2 

= (7-2) 

k ct=1 
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If we restrict the field to the interior of a 1-dimensional (but large) volume V then the normal modes of the field are 
simply 

u\{x) = (2Vuj k )- 1/2 coskx (7.3) 
u 2 k (x) = (2Vuj k )- 1/2 sinkx (7.4) 

so that the time dependence of the fc th mode has the following representation in terms of creation and annihilation 
operators 

^(t)=a fc e- iu ** + 4e*-* t (7.5) 
<p 2 k (t)=i(a k e- i "* t -ate i "* t ) (7.6) 
In terms of this mode decomposition, the Hamiltonian is 

H = \ E E («) 2 + e ^ 2 ) + ^n + \ mn2q2 - A E E «(*K( ) ^ (*)■ ( 7J ) 

ct— 1 a— 1 

Notice that the coupling constant A in the last term can be grouped with the mode function u£(0) to give an effective 
coupling constant that depends on the particular mode C£ = Att^(O). Therefore, the Hamiltonian for this system is 

H = \ E E («) 2 + fc2 (^) 2 ) + £ + ^ m "V -EE CfcffWvsw = ^-Ho qbm, (7.8) 

fc <t=1 fc <r=l 

which is precisely the Hamiltonian for a harmonic oscillator q(t) coupled to an environment composed of a bath 
of harmonic oscillators In other words, the MOF model for a mirror at rest can be related to quantum 

Brownian motion where the field provides the environment that the mirosc interacts with. QBM has a long history 
and is well-studied so that results already found in that literature can be applied directly to the interaction of a field 
with a static mirror via the MOF model. For example, the master equation is exactly known for this system [68] and 
so one can study its behavior near the perfectly-reflecting limit where A — > oo or, equivalently, m — > as well as in a 
non-zero temperature regime. 

A similar result holds for TV mirrors held at rest at positions x = L a with a = 1, . . . , N. It is straightforward to see 
that the corresponding Hamiltonian, when decomposing the field into harmonic oscillators, is 

1 2 N / 2 1 2 \ 

H = 2 E E («) 2 + fc2 (^) 2 ) + E + 2 TO ^ 2 "EE C kaVa(tM(t) = ffiV-HO QBM (7.9) 

h a=\ o=l \ a k cr=l ' 

where the effective bilinear coupling constant is C£ a = Au%(L a ). Therefore, N static mirrors in the MOF model 
correspond to N harmonic oscillators (mirosc variables) coupled to a bath of oscillators (the field). For N = 2 
oscillators coupled to a general environment, the exact master equation has been derived in [69] and thus can be used 
to provide a different perspective and new insights in the description of a field coupled to two partially transmitting 
mirrors via the MOF model. 

B. Slowly moving mirrors and QBM 

Turn next to find the relationship between slowly moving mirrors in the MOF model and quantum Brownian 
motion. Let us first consider one mirror since the result for N mirrors will generalize in an obvious way. Assume that 



the mirror is in an externally generated potential V(x), such as a harmonic trap. Then the Hamiltonian in (7.7) is 

„ h -mWq 2 + — — 

2m 2 H 2M 



H = \ E E («) 2 + fc2 (^) 2 ) + £ + l mn2 i 2 + vrr + V ^ A E E <K*K(s(*)) m) (7.io) 



where we have included the worldline variable to the Hamiltonian. Notice that from a QBM perspective, the effective 
coupling constant acquires a time dependence since the mode function is now time dependent, u%(Z(t)). However, if 
the potential V(x) restricts the motion of the mirror to be only small perturbations from its equilibrium position at 
x = then we may expand the mode function about the origin so that the interaction term above becomes 



2 



A E E 9(t)u%{Z(t)) <p%(t) = - £ E C%q(t)<p%(t) - \Z(t) £ £ d x u%(0)q(tM(t) + 0(Z 2 ) (7.11) 



k 



23 



Therefore, the Hamiltonian (7.10) is equal to an unperturbed Hamiltonian, given by the 1-harmonic oscillator QBM 



Hamiltonian in (7.8 1, plus an interaction Hamiltonian that describes perturbations due to the small displacement of 



the mirror that arise from interactions between the field oscillators and the mirosc, 

2 

H = Hi-ho qbm - XZ(t) ]T d x u%(0)q(tM(t) + 0{Z 2 ) (7.12) 



k ff=l 



Hence, one can compute the perturbations of, for example, the exact master equation for 1-harmonic oscillator QBM 
to study the behavior of a movable, partially transmitting mirror. Notice that if V(Z) = A/il§Z 2 (i)/2 then (7.12) 



describes a nonlinearly coupled QBM system where the mirosc and the mirror's position are the two oscillators in an 
open system that couples to the bath provided by the field oscillators. The nonlinearity is only in the mirror's position 
(i.e., from the 0{Z 2 ) terms above) but the mirosc and the field oscillators still couple to each other bilinearly. 
The generalization to N mirrors should be obvious with the Hamiltonian describing the system being 



N 



H = H N -no qbm - A]T Z a (t) E d x ul(L a )q(tM(t) + 0{Zl) (7.13) 

a— 1 k (7—1 

where the unperturbed position of the a th mirror is at x = L a . In particular, one can compute the perturbations of, 
for example, the exact master equation for 2-harmonic oscillator QBM to study entanglement, decoherence, etc., 
of a cavity with movable, partially transmitting mirrors. 



VIII. SUMMARY AND FURTHER DEVELOPMENTS 



In this paper we constructed a microphysics model of moving mirrors interacting with a quantum field. The novel 
ingredient we introduced is a harmonic oscillator (a "mirosc") model describing the internal degrees of freedom of 
the mirror that couples to the incident radiation thereby providing a mechanism for the dynamical interplay of the 
mirror-field system. Since the field can transfer (receive) energy and momentum to (from) the mirosc the collection 
of them serves the function of a partially reflecting or transmitting mirror. We showed that this mirror-oscillator-field 
(MOF) system can perfectly reflect or perfectly transmit radiation depending on the values of the mirosc mass m, 
natural frequency f2, and coupling strength A to the field. Perfect reflection can be attained in three ways: 1) m — > 0; 
2) A — > oo; and 3) Frequency uj of an incident wave is equal to the mirosc natural frequency J7. Limits 1) and 2) 
exhibit perfect reflection (or nearly so) among a broad frequency bandwidth whereas limit 3) strongly reflects modes 
with frequencies near !7 because of a resonant excitation of the mirosc. 

The MOF model reduces to several commonly used models of moving mirrors in a quantum field. We showed that 
when the mirosc variable q(t) evolves adiabatically (|<jr| <C \Q 2 q\) or when m — > but mil 2 = k remains constant then 
the MOF model reduces to the Barton and Calogeracos (BC) model [2] of a partially transmitting moving mirror. The 
free parameter in the BC model 7 is related to the mirosc parameters of the MOF model (m, fi, A) by 7 = A 2 /(2mil 2 ). 
The "auxiliary field" model of Golestanian and Kardar [151135] arises from the MOF model in the limit that m — » 0. In 
this limit, there is no mirosc dynamics and q(t) becomes an auxiliary variable. In the quantum theory, q(t) may have 
any possible realization (see (2.29)), which manifests as a Dirichlet boundary condition on the field at the location of 
the mirror and thus perfectly reflects incident radiation (see (2.31)). We also showed that our MOF model reduces 
to the phenomenological model of a mirror interacting with a cavity field via the radiation pressure exerted on the 
mirror's surface when a number of assumptions are made (though these may not all be necessary to derive the Nx- 
coupling in other setups) . This "TVcc-coupling" is often used to describe laboratory setups but may be extended by the 
MOF model to scenarios where the mirosc does not evolve adiabatically, which may exhibit interesting macroscopic 
(or perhaps mesoscopic) quantum phenomena. Additionally, iVa;-type coupling provides the leading order corrections 
to the classical radiation pressure coupling when the cavity is occupied by low photon numbers. The model we present 
in this paper will remain useful even when the necessary conditions for it to match with models with iVx-type coupling 
are not met, for example, when the mirror motion is sufficient to excite field quant to higher modes. 

The bulk motion of the mirror in the MOF model, which may be relativistic depending on the application, can be 
derived from an action or a Hamiltonian. In either formulation, we find that the mirror moves with a time-dependent 
effective mass M e g that is composed of the mirror's rest mass M and the mirror's total internal energy, which comes 
from the energy of the mirosc itself and its interaction with the field. We also showed (in a purely classical setting) 
that the MOF model seems to admit physical solutions despite the use of a point particle description for the mirror's 
motion and despite the interaction between the mirror and field resembling that of a charged particle (which can be 
plagued by pathologies). We demonstrated that when the field is generated by its interaction with the mirosc alone 
so that there is no external source J cxt and no initial field configuration present then the mirror will evolve on an 
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inertial trajectory (i.e., constant velocity), which is the correct expected result, in contradistinction to the radiation 
reaction on a point charge in electrodynamics where the charge may exhibit run-away motions in the absence of any 
external forces acting on the charge. 

As an application of the MOF model, we studied the "cooling" of a mirror by its interactions with an external field 
in a purely classical context. We found that when the mirosc is weakly coupled to the field that the mirror, when 
perturbed, will oscillate around its equilibrium configuration while its displacement amplitude decays slowly in time. 

An interesting consequence of our MOF model for moving mirrors is that it relates to models of quantum Brownian 
motion (QBM) in a straightforward manner. The relation essentially follows because the field can be regarded as a 
continuum of harmonic oscillators. Hence, for N mirrors held at rest, the MOF model is equivalent to N harmonic 
oscillators in a bath of oscillators (from the field). For N — 1,2, the master equation for such a system in a general 
environment has been derived exactly [551 ISS] an d even for general N [70] . Consequently, the MOF model can be used 
to study the superposition of two mirrors, the decoherence by and the disentanglement of moving mirrors via a field, 
etc., so as to gain insight into these aspects of macroscopic quantum phenomena. We expect that the rich repository 
of technical tools and physical insights from the study of QBM can be carried over directly to our MOF model for a 
broad range of applications involving moving mirrors and quantum fields. For example, QBM results for systems at 
finite temperature may provide a simple way to incorporate thermal effects into the MOF model. We will begin to 
explore this theme in a follow-up paper [53] on the theory of OM from an open quantum system viewpoint. 

The generalization of the MOF model to 3 spatial dimensions can be made where the mirror is an extended body 
having some surface geometry. On this surface, we may place a layer of mirosc's that play the role of the electrons in a 
metal gas or dielectric medium providing the mirror's light degrees of freedom and responsible for reflection of incident 
radiation over some bandwidth of the electromagnetic spectrum (e.g., optical as in many metals). Incorporating the 
electromagnetic field in the MOF model should also be straightforward as its structure is similar to that of a minimally- 
coupled scalar field in the MOF model (see the appendix of [61]). 

In the second series on back-action effects we will study the full quantum mechanical evolution of the MOF system 
in the context of mirror cooling. Therein, we will derive the exact equations of motion describing the mirror's average 
position. In the most general case we will show that the mirror motion is described by an integro-diffcrcntial equation 
exhibiting non-Markovian dynamics. The equations can be simplified through a series of approximations which directly 
relate to experimentally engineerable quantities, such as the cavity's quality factor, and the relevant timescales for the 
mirror's internal dynamics. Given the broad range of applicability, these results can be employed to guide theoretical 
and experimental investigations ranging from the cooling of the center of mass motion of moveable mirrors, having 
broad-band reflective properties, to the manipulation of trapped ions near surfaces, possessing narrow-band reflective 
properties |14j . 




In the third series we will address the moving mirror analog of the back-reaction of Hawking radiation [19] on 
the evolution of a black hole. There are controversies in some deep issues related to the end-state of black hole 
evaporation resulting from the Hawking effect, namely, whether complete evaporation of a black hole means the non- 
unitary evolution of quantum states (see, e.g., [71] ) which violates the basic tenets of theoretical physics or if unitarity 
is preserved, and if so, how? One key ingredient, the back-reaction of the emitted radiation on the spacetime, has 
not been taken into account fully or correctly (for a recent update, see [72] [73] and papers cited therein.) There are 
analog studies on how information is shared in the black hole (harmonic) atom - quantum field system (see, e.g., [74] 
and references therein.) as well as moving mirror analog problem |17| . The connection was made between the s-wave 
component of Hawking evaporation and the emission of radiation from moving mirrors by the dynamical Casimir 
effect but, like the original calculation by Hawking, treated the effects of back reaction rather coarsely. Since the 
MOF model offers a large degree of flexibility and tractability, we were able to find exact equations of motion for 
the mirror incorporating the effects of back- reaction [BUJ. These exact solutions, as well as those from the atom-field 
analogs, can provide new insights into this basic issue in theoretical physics. 
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